Trapped-Ion Quantum Simulator: Experimental Application to Nonlinear 

Interferometers 
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We show how an experimentally realized set of operations on a single trapped ion is sufficient to 
simulate a wide class of Hamiltonians of a spin-f/2 particle in an external potential. This system 
is also able to simulate other physical dynamics. As a demonstration, we simulate the action of an 
n-th order nonlinear optical beamsplitter. Two of these beamsplitters can be used to construct an 
interferometer sensitive to phase shifts in one of the interferometer beam paths. The sensitivity in 
determining these phase shifts increases linearly with n, and the simulation demonstrates that the 
use of nonlinear beamsplitters (n=2,3) enhances this sensitivity compared to the standard quantum 
limit imposed by a linear beamsplitter (n=l). 



One of the motivations behind Feynman's proposal 
for a quantum computer [1] was the possibility that one 
quantum system could efficiently simulate the behavior of 
other quantum systems. This idea was verified by Lloyd 
[2] and further explored by Lloyd and Braunstein [3] for a 
conjugate pair of variables such as position and momen- 
tum of a quantum particle. Following this suggestion we 
show below that coherent manipulation of the quantized 
motional and internal states of a single trapped ion us- 
ing laser pulses can simulate the more general quantum 
dynamics of a single spin-1/2 particle in an arbitrary ex- 
ternal potential. Previously, harmonic and anharmonic 
oscillators have been simulated in NMR [4] . 

In addition to demonstrating the basic building blocks 
for simulating such arbitrary dynamics, we experimen- 
tally simulated the action of optical Mach-Zehnder inter- 
ferometers with linear and nonlinear second- and third- 
order beam-splitters on number-states. Interferometers 
with linear beamsplitters and nonclassical input states 
have engendered considerable interest, since their noise 
limits for phase estimation can lie below the standard 
quantum limit for linear interferometers with coherent in- 
put modes [5-8] as has been demonstrated in experiments 
[9]. A number of optics experiments have exploited the 
second-order process of spontaneous parametric down- 
conversion [10], which can be regarded as a nonlinear 
beamsplitter. By cascading this process, a fourth-order 
interaction has also recently been realized [11]. One diffi- 
culty in these experiments is the exponential decrease in 
efficiency as the order increases, necessitating data post- 
selection and long integration times. In the simulations 
reported here, nonlinear interactions were implemented 
with high efficiency, eliminating the need for data post- 
selection and thereby requiring relatively short integra- 
tion times. 

To realize a quantum computer for simulating a spin 
s = 1/2 particle of mass ^ in an arbitrary potential, one 



must be able to prepare an arbitrary input state 
[^'(m,,^)) =5Z(ci-l i)\n)+c^n\ T)|ri)), 



(1) 



where the particle's position wavefunction is expanded 
in energy eigenstates |n) of a suitable harmonic oscillator 
and \ms) [nis £ {[, t}) represent the spin eigenstates in a 
suitable basis. We have recently demonstrated a method 
to generate arbitrary states of the type in Eq.(l) in an 
ion trap [12,13]. The computer should then evolve the 
state according to an arbitrary Hamiltonian 
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where we require only that the potential V{z, rris) can be 
expanded as a power series in the harmonic oscillator lad- 
der operators a and a) and be approximated to arbitrary 
precision by a finite number of terms with maximum or- 
der TV. The nis are a set of observables in a general 
two-level Hilbert space that can all be mapped to a lin- 
ear combination of the indentity / and the Pauli matrices 
Uj . The operators a± are defined as a± = i iuy , all 
/3„m,X„m are complex numbers, and all a„m,7„m are 
real numbers. 

In our realization of an analog quantum computer we 
consider the Hamiltonian of a trapped atom of mass /z, 
harmonically bound with a trap frequency uJz and inter- 
acting with two running-wave light fields having a fre- 
quency difference Aw and a phase difference ip at the po- 
sition of the ion. Both light fields are assumed to be de- 
tuned from an excited electronic level so they can induce 
stimulated-Raman transitions between combinations of 



1 



two long-lived internal electronic ground-state levels with 
energy difference tiujo and the external motional levels of 
the ion [14]. For our purpose it is sufficient to consider 
the motion along one axis in the trap. After applying a 
rotating wave approximation and adiabatic elimination 
of the near resonant excited state [14] , and switching to 
an interaction picture of the ion's motion, the resonant 
interaction for Raman beam detuning Alv = ewo + I'jJz 
(e = {0, 1}, I integer) can be written in the Lamb-Dicke 
limit {r]'^{{a + a^f) <C 1) as [15,16] 



have not been previously demonstrated. One of the ex- 
periments discussed below uses two ffis pulses, therefore 
demonstrating the feasibility of generating Hq^ as well 
[21]. 

As a demonstration of quantum simulation using a sin- 
gle trapped atom, we employ the interactions Hu,Hi2, 

and His to efficiently simulate a certain class of n-th or- 
der optical beamsplitters described by Hamiltonians 



^l,\l\^JTJ- + - %\l\) 1,1, 



h.c 



(3) 



The coupling strength Q, is assumed to be small enough 
to resonantly excite only the l-th spectral component. 
The Lamb-Dicke parameter r] = Ak zq is the product of 

the z-projection of the wavovcctor-diffcrcncc Afc of the 
two light fields and the spatial extent of the ground state 
wave function zq = \J h, j i^mu: ^ . For e = 1 the internal 
state changes during the stimulated Raman transition 
and the interaction couples | i)|n) ^ | t)|'^ + 0; while 
for e = only motional states |n) <-^- \n + l) are coupled 
with a strength independent of the internal state [17]. 

Following Lloyd and Braunstein [3,18], by nesting and 
concatenating sequences of H^i operations according to 
the relation 

e-*^**e-*^'**e*^**e*^'** = gF^I^'^'l**' + 0{bt\ (4) 

the set of operators {i?oi, -^02, -f^oa, -f^io, -^ii, ^^12, ^13} 
is sufficient to efficiently generate arbitrary Hamiltoni- 
ans. This conclusion is straightforward for the spin, since 
{cr+, CTz} are a complete basis of that algebra. For in- 
teractions that only involve the motion (e = 0) it follows 
from the fact that 

[i?02, -ffos] i \a.(^c? -\- a*(a^)^a} -|- lower orders (5) 



and 



aa' a 



t«2 



a*(at)2a,/3(a^)"a" + /3*(a^)'"a"] = 



(2m - n) [a/3(a'f)'"a"+i a/3*(at)"a'"+i - /i.e.] 

-h lower orders, (6) 

so one can build up arbitrary orders in the effective 
Hamiltonian by recursive use of Eq.(4). Similar argu- 
ments hold for the set of interactions, and by com- 
bining both types of interactions, the series expansion of 
the Hamiltonian in Eq.(2) can eventually be constructed. 

Most of these interactions have been demonstrated in 
previous ion-trap experiments. i?io is usually called the 
carrier interaction, i/oi and i7o2 arc coherent and squeeze 
drives respectively and H\\,H\i are first and second blue 
sideband [19,20]. The third-order interactions iJosi-f^is 



B„ = ?if^„[a(6t)" + at(6)n]. 



(7) 



Here a and h are the usual harmonic oscillator lowering 
operators for the two quantized light modes, r2„ is the 
coupling strength, and we simulate the special case where 
the number of photons in mode a is or 1 and n = 1 , 2 
or 3. Two such beamsplitters can be used to construct 
a Mach-Zehnder interferometer as sketched in Fig. 1. 
The order n = \ corresponds to the commonly used lin- 
ear beamsplitter that is typically realized by a partially 
transparent mirror in experiments. Such interferometers 
can measure the relative phase of the two paths of the 
light fields that are split on the first beamsplitter and 
recombined on the second. The phase can be varied by 
changing a phase shifting element (the box labeled </) in 
Fig. 1) and detected (modulo 27r) by observing the in- 
terference fringes of the recombined fields. We restrict 
our attention to a pure number-state \n = l)a imping- 
ing on the first beamsplitter from mode a and a vacuum 
state |n = 0)b from mode b. After propagating the input 
state through the first beamsplitter with f7„ adjusted to 
give equal amplitude along the two paths in the output 
superposition, the state becomes 

|l>a|0>6-^(|l>a'|0)6' + |0)„Hn>6')- («) 

Phase shifters in optical interferometers alter a classical- 
like coherent state \a) to one that is shifted to \ae^'^). 
In the context of Fig.l this phase shift corresponds to 
\n) — > e'"'^|n) for a number state, leading to 
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-J={\l)a-P)b' + P)a'\n)y)^ 

-^(|l)„HO)6' + e-10V|n)60. 



(9) 



The second beamsplitter recombines the two field modes 
leading to an average probability of 



("a") = ^[l-cos(n(^)] 



(10) 



for detecting one photon in the output arm with the de- 
tector in Fig.l. 

We have experimentally simulated the nonlinear beam- 
splitter of Eq.(7) using a single trapped ^Be"*" ion. The 
operator a is replaced by cr"*", the raising operator be- 
tween two hyperfine states |F = 2,mF = —2) = |l)a 
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and \F = l,niF = —1) = \0)a in the ground state 

manifold. These operators (and also their respective Her- 
mitian conjugates) are not strictly equivalent, but their 
action is the same as long as we restrict our attention 
to situations that never leave the {|0)a,|l)a} subspace. 
The simulated linear and nonlinear interferometers fulfill 
this restriction, as long as the input state is |l)a|0)5. The 
optical mode with lowering operator b is replaced by the 
equivalent harmonic oscillator mode of motion along one 
axis in the trap, with number states \n). 

Our experimental system has been described in de- 
tail elsewhere [19,20,22]. We trapped a single ^Be+ ion 
in a linear trap [23] with motional frequency = 2tt 
3.63 MHz (Lamb-Dicke parameter rj =0.35) and cooled 
it to the ground state of motion. The trap had a heating 
rate of 1 quantum per 6 ms [23] that was a small pertur- 
bation for the duration of our experiments (< 260 /xs). 
After cooling, the ion was prepared in the |l)o|0)6 state 
by optical pumping. Starting from this state we used 
Raman-transitions to drive a 7r/2-pulse on the ion's n-th 
blue sideband {Hj cx (t+(5^)" + h.c), creating the state 
( 1 1) a' 1 0) fc' + 1 0) a' I n) 6/ ) / V2. For different orders n the 7r/2- 
pulse time scales as Vn^-frj" [14]. The observed 7r/2-times 
of (4.0, 17.3, 115) fis do not scale exactly as the theoretical 
prediction due to different laser intensities used for the 
different values of n. A phase shift 4> = Awz t was then 
introduced by switching the potential of the trap endcaps 
to a different value for time i, thus changing the motional 
frequency by a fixed amount Auz- After a second 7r/2- 
pulse on the n-th sideband we measured the probability 
{ua") for the ion to be in |l)a- The interference fringes 
created by sweeping t are shown in Fig. 2. The final 
state of the ion oscillated approximately between \l)a" 
and |0)a" as t was varied, with frequency nAujz- 

In interferometric measurements, we want to maxi- 
mize our sensitivity to changes of (j) around some nominal 
value. We therefore want to minimize 

S4>= (11) 

\d{na")/d(j)\ 



where AA = y {A^) — {A)"^ is a measure of the fluctua- 
tions between measurements of an operator A. Eq. (11) 
applies to our simulator with (j) = AtUzt. In our experi- 
ments 

(na") = (C/2) [l-cos(n0)], (12) 

where C is the contrast of the observed fringes. Ideally 
C = 1 [Eq. (10)] but is observed to be < 1 due to im- 
perfect state preparation and detection, and fluctuations 
in the ambient magnetic fleld and the trap frequency . 
The sensitivity of the interferometer is max;imized when 
the slope of {fia") with respect to (j), d{ha") / d4> is max- 
imized, that is, for values of where ncj) = ■7rfc/2, k an 
odd integer. We characterize the fluctuations Aua" with 



the two-sample Allan variance, commonly used to char- 
acterize frequency stability [24]. In the present context, 
a series of M (total) measurements oiha" = |l)o"(l|a" is 
divided into bins of Nb measurements averaged according 
to 

ii+l)Nt-l 

{ha")i = l/Nb (13) 

where 2 < Ni, < M/2 and {na")j is the j-th measurement 
of fia"- The Allan variance characterizing fluctuations 
between measurements is given by 

{(Tn^„{Nb)f = _ {{na")i+l - {na")if- 

(14) 

Making the identification an „ (A"b) = Ana", in Fig. 3, 
we plot 6(j) vs. iVf,. The solid curve is the theoretical 
standard quantum limit for a linear interferometer with 
perfect contrast and unity detection efficiency, given by 
Aua" /\/Nb where (Ana")^ is the variance due to projec- 
tion noise [25]; Ana" = 0.5 at the points of maximum 
slope in our fringes. The simulation of the linear inter- 
ferometer shows only a small amount of excess noise over 
the theoretical limit, due mainly to the C = 0.92 con- 
trast of the fringes, while the nonlinear interferometer 
simulations have a noise-to-signal ratio below the linear 
interferometer standard limit. The potential gain in slope 
for n = 3 is almost exactly cancelled by the loss in fringe 
contrast, so the noise-to-signal ratio for n = 2 and n = 3 
is about the same. 

In conclusion, we have shown how coherent stimulated- 
Raman transitions on a single trapped atom can be used 
to simulate a wide class of Hamiltonians of a spin- 1/2 par- 
ticle in an arbitrary external potential. This system can 
also be used to simulate other physical dynamics. As a 
demonstration, we have experimentally simulated the be- 
havior of n-th order nonlinear optical beam-splitters act- 
ing in a restricted Hilbert space. Our simulation demon- 
strates how interferometer sensitivity improves with the 
order of the beam splitter. As a practical matter, the 
2nd- and 3rd-order beamsplitters demonstrated here give 
increased sensitivity for diagnosing motional frequency 
fiuctuations in the trapped-ion system. With anticipated 
improvements in motional state coherence [23], it should 
be possible to simulate more complicated Hamiltonians. 
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FIG. 1. Principle of the Mach-Zehnder interferometer. 
Modes 'i! a and are superposed on the first beamsplitter. 
After the beamsplitter has acted the modes a' , ^6' are prop- 
agated along separate paths to a second beamsplitter. Mode 
^/f,/ may undergo a variable phase shift induced by the phase 
shifter 0. Modes 'i! a" and ^yi emerge after the second beam- 
sphtter and one of the modes is put onto a detector. Varying 
(j) will lead to a sinusoidal behavior of the intensity on the 
detector (fringes). 




(a) 

0=0.92 



FIG. 2. Interference fringes for simulated interferometers 
(a) of order n = 1, (integration time per data point was 0.50 
s); (b) n = 2, (0.53 s) and (c) n = 3, (0.63 s). {ria") is 
the probability to find the ion in |l)a", while t is the time 
for which the trap frequency was shifted by AtJz, directly 
proportional to the phase shift cj> — Aujz t. The frequency of 
the fringes increases linearly with order n. C is the observed 
contrast of the fringes. 
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FIG. 3. Noise-to-signal ratios 5(j){Nb) for the n — 1 linear 
interferometer (solid squares), n = 2 nonlinear interferome- 
ter (solid circles) and n = 3 nonlinear interferometer (open 
circles) vs. the number of measurements A'^b. The solid line 
is the theoretical limit for the linear (n=l) interferometer, 
assuming perfect contrast and detection efficiency. 



